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摘要 : ”从 RNSS (radio navigation satellite system) 米 级 服务 角度 给 出 了 BDCS (Beidou coor- 
dinate system) 实现 精度 估算 方法 ， 包 括 广播 星 历 定 位 法 、 精 密 星 历 与 广播 星 历 比 对 法 和 卫星 激 
光 测 距 (satellite laser ranging, SLR) 坐标 估计 法 。 广 播 星 历 定位 法 适用 于 用 户 对 参考 框架 实现 
= 精度 的 实时 监测 ， 但 监测 精度 受 限 于 广播 星 历 精度 。 轨 道 比 对 法 引入 了 事后 精密 产品 ， 监 测 精 度 
> 较 高 ， 但 其 时 效 性 相对 较 差 。SLR 坐标 估计 法 不 受 电离 层 延迟 和 卫星 钟 差 影响 ， 在 三 种 方法 中 
精度 最 高 ， 但 由 于 激光 测 距 数据 量 较 少 ， 无 法 做 短期 监测 。 结 果 表 明 ， 三 种 方法 估算 BDOS 与 
ITRF 对 齐 精度 均 为 厘米 级 。 最 后 ， 在 中 俄 卫 星 导航 系统 协作 背景 下 ， 估 算 了 BDCS 和 PZ-90 
对 齐 精度 ， 结 果 表 明 ， 两 者 在 厘米 量 级 对 齐 ， 能 够 满足 米 级 定位 用 户 应 用 需求 。 随 着 北斗 三 号 / 
播 星 历 精 度 的 提升 和 SLR 的 发 展 ， 上 述 三 种 方法 的 估算 精度 将 会 进一步 提高 。 
关 键 词 : 北斗 坐标 系 ; 广播 星 历 定位 ;卫星 激光 测 距 ， 精 度 监测 
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1 引 言 


fez 


国际 地 球 参考 框架 (international terrestrial reference frame, ITRF) 是 国际 地 球 参考 系 
(international terrestrial reference system, ITRS) 的 实现 ， 通 过 全 球 导 航 卫 星系 统 (global 
= navigation satellite system，GINSS)、 甚 长 基线 干涉 测量 (very long baseline interferometry, 
FT VLBI) SLR 以 及 星 载 多 普 勒 定 轨 和 无 线 电 定位 (doppler orbitography and radio-positioning 


integrated by satellite, DORIS) 四 种 空间 大 地 测量 技术 建立 和 维持 。2016 年 1 月 ITRF 发 

布 最 新 版 国际 地 球 参考 框架 的 ITRF2014 “。 
不 同 的 导航 系统 采用 不 同 的 空间 基准 。 如 GPS 采用 WGS84 坐标 系 ，GLONASS 采 

PZ-90 坐标 系 ，Galileo 采用 GTRF 坐标 系 。 为 保证 各 卫星 导航 系统 间 的 兼容 和 互 操作 ， 


前 各 导航 系统 参考 框架 采用 向 ITRF 对 齐 的 方式 实现 导航 系统 之 间 空 间 基准 的 一 致 性 。 对 
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航 系统 参考 框架 的 精度 可 以 采用 两 种 方法 进行 估算 : (1) 利用 地 面 站 在 两 个 参考 框架 的 精密 
坐标 之 间 直 接 计算 转换 参数 ;(2) 通过 包含 坐标 框架 信息 的 卫星 轨道 产品 和 定位 结果 间接 监 
测 估算 框架 对 齐 精度 。 对 于 导航 系统 供应 商 ， 一 般 采 用 第 一 种 方法 获得 与 ITREF 对 齐 的 精确 
估算 。 该 方法 需要 参考 框架 实现 过 程 中 选用 的 监测 站 数据 ， 而 该 数据 一 般 用 户 不 可 获得 。 为 
实现 对 导航 系统 参考 框架 精度 的 监测 ， 一 般 采 用 第 二 种 方式 ， 即 通过 广播 星 历 等 产品 间接 估 
算 。 两 种 方法 的 监测 精度 和 时 效 性 有 所 区 别 : 前 一 种 方法 可 实现 毫米 级 精度 估算 ， 一 般 用 于 
框架 实现 时 的 估算 ; 第 二 种 方法 精度 较 第 一 种 低 ， 但 更 方便 用 户 对 参考 框架 精度 进行 监测 和 
估算 。 本 文 重点 分 析 第 二 种 方法 的 估算 精度 ， 并 利用 实测 数据 对 BDCS 开展 评估 ”。 

北斗 卫星 导航 系统 的 空间 基准 是 北斗 坐标 系 。 北 斗 坐标 系 采 用 120 余 个 分 布 于 中 国境 内 
及 全 球 的 GNSS 监测 站 观测 数据 并 实现 与 ITRF2014 的 对 齐 。 最 新 版 北斗 坐标 系 通过 IGS 
监测 站 实现 与 ITRF2014 对 齐 ， 对 齐 精度 约 为 毫米 级 。 

除 导 航 系 统 公 布 的 导航 系统 参考 框架 的 精度 ，IGS 从 用 户 角度 也 提供 GPS 空间 基准 的 
精度 监测 结果 。IGS 通过 比较 GPS 广播 星 历 与 IGS 快速 星 历 产品 比 从 而 计算 7 参数， 实现 
对 空间 基准 对 齐 精 度 监 测 ， 在 一 定 程度 上 代表 了 GPS 卫星 对 RNSS 用 户 的 空间 基准 传递 精 
pee = 

E E E tee E 
及 所 有 卫星 搭载 激光 角 反 射 器 的 特点 ， 广播 星 历 定位 方法 以 及 SLR 坐标 方法 开展 对 
BDS 空间 基准 传递 精度 进行 估算 。 eee sx 间 基 准 互 操作 应 用 背景 ， 利 用 SLR 方法 估算 
T BDS 5 GLONASS 系统 空间 基准 对 齐 精 度 ， 为 导航 系统 兼容 互 操作 应 用 提供 支撑 。 


2 参考 框架 实现 精度 监测 方法 


2.1 ”广播 星 历 定位 法 
广播 星 历 定位 法 的 处 理 采 用 固定 广播 星 历 卫 星 轨 道 和 钟 差 ， 估 算 GNSS 监测 站 每 天 一 
组 位 置 参数 、 接 收 机 各 历 元 钟 差 、 每 6 h 大 气 折射 误差 和 相位 模糊 度 。 式 (1) 为 观测 方程 : 


(1) 
LC; = p+ c(t, | ts) | dtrop Piet Ni + dother ， 


{re =p ote c(t, + ts) lr dtrop T dother ’ 


AF, PO; 为 伪 距 无 电离 层 组 合 值 ，LC; 为 相位 无 电离 层 组 合 值 ，p 为 星 地 间 几 何 距离 ; tr 
All ts 分 别 为 接收 机 钟 差 和 卫星 钟 差 ，duop 为 对 流 层 延 迟 ，dotne 为 其 他 误差 项 ， 如 相对 论 效 
应 ， 国 体 潮 海潮 改正 以 及 多 径 噪 声 等 ， 和 Xi 为 无 电离 层 组合 频 率 载波 波长 ，N; 为 组 合 模糊 度 
未 知 数 ;，c 为 光速 ，c = 299792 458 m .s-1。 为 达到 对 导航 系统 空间 基准 传递 实时 监测 估算 
的 目的 ， 广 播 星 历 定位 法 采用 广播 星 历 计算 各 历 元 卫星 轨道 和 钟 差 。 

利用 IGS 台 站 一 天 弧 长 观测 数据 估算 坐标 参数 ， 解 算出 监测 站 在 导航 系统 参考 框架 下 
的 坐标 。 利 用 上 述 坐 标 及 该 监测 站 在 ITRF2014 框架 下 的 坐标 估算 转换 参数 ， 即 得 到 该 参考 
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框架 相对 于 ITRF2014 的 平移 参数 、 旋 转 参 数 和 尺度 因子 。7 参数 转换 公式 如 下 : 


XaNss = X trans + (1 + Xscale) X X rot X XITRF ) (2) 


式 中 ，Xuans 为 平移 参数 ，Xiot 为 旋转 参数 ，X,a 为 尺度 因子 ，Xess 为 广播 星 历 定位 法 
得 到 的 监测 站 坐标 估计 值 ，Xirar 为 监测 站 在 ITRF2014 框架 下 的 坐标 。 
2.2 SALEM 

广播 星 历 空间 基准 为 各 导航 系统 参考 框架 ， 精 密 星 历 空间 基准 主要 由 IGS 维持 。 目 前 
IGS 空间 基准 与 ITRF2014 一 致 ， 因 此 可 以 通过 广播 星 历 计算 导航 卫星 的 位 置 与 精密 星 历 进 
行 比较 ， 获 得 导航 系统 参考 框架 与 ITRE 的 转换 关系 。 由 于 IGS 的 精密 轨道 是 卫星 质心 坐 
标 ， 而 广播 星 历 为 卫星 天 线 相位 中 心 位 置 ， 因 此 需要 进行 卫星 天 线 相 位 中 心 到 质心 的 改正 。 
考虑 到 北斗 不 同类 型 卫星 姿态 控制 模式 不 同 ， 卫 星 天 线 相 位 中 心 修正 具体 如 下 。 

(1) IGSO/MEO 卫星 天 线 相 位 中 心 改正 模型 ; 


r Tsun — T 
; E a ez = Cy X ez ; (3) 
sun 


HP, r, rsun 为 惯性 系 中 卫星 位 置 向 量 和 太阳 位 置 向 量 ， ej, €y, ez 为 星 固 系 到 地 固 系 的 转 
换 和 矩阵。 
(2) GEO 卫星 天 线 相位 中 心 改正 模型 : 


VU 
ez 一 一 一 ，ey 一 ez X ，en 一 ey Xe , (4) 
Ir] E |v] i 
Ht, r, v 为 惯性 系 中 卫星 位 置 向 量 和 速度 向 量 。 
相位 中 心 改正 模型 为 : 
X Zphs 
Y| = Rciscts x lex €y ez X | Yphs ) (5) 
Z Zphs 


EH, Resets 为 惯性 系 与 地 固 系 转换 旋转 个 阵 ，zpps, Ypns> Zpns 为 星 固 坐标 系 下 卫星 天 线 相 
位 中 心 偏差 ，X, Y, Z 为 卫星 天 线 相位 中 心 改正 值 。 

归 算 到 同一 中 心 后 ， 计 算 互 差 并 进行 7 参数 转换 ， 即 得 到 该 参考 框架 对 齐 ITRF 的 平移 
参数 、 旋 转 参数 和 尺度 因子 。7 参数 转换 公式 同 公式 (2)。 
2.3 SLR 坐标 估算 法 

SLR 坐标 估算 法 与 广播 星 历 定位 法 类 似 ， 通 过 将 SLR 台 站 在 导航 系统 参考 框架 中 的 坐 
标 与 ITRF 坐标 进行 比较 ， 获 得 两 个 参考 框架 的 转换 参数 。 激 光 数 据 为 地 面 站 到 卫星 角 反 射 
器 之 间 的 测 距 ， 具 有 受 电离 层 折射 误差 影响 小 、 测 量 精度 高 、 且 不 受 钟 差 影响 等 特点 。 北 斗 
三 号 卫星 基于 星 间 链 路 技术 提供 实时 全 球 高 精度 广播 星 历 产品 ， 轨 道 精度 高 "。 同 时 ， 北 斗 
导航 系统 各 卫星 均 配 备 了 角 反 射 器 ， 因 此 可 以 采用 广播 星 历 轨道 产品 和 SLR 高 精度 测量 数 
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据 获得 SLR 台 站 在 BDCS 框架 下 的 坐标 。 将 该 坐标 与 SLR AMEE ITRE 框架 下 的 坐标 进行 
比较 ， 可 获得 两 参考 框架 的 转换 参数 ””。 

将 广播 星 历 计算 的 轨道 由 相 心 坐标 改正 至 质心 坐标 ， 方 法 同上 ， 并 进行 卫星 质心 改正 。 
采用 一 个 月 观测 弧 长 对 SLR 台 站 坐标 进行 估算 ， 并 将 估算 坐标 与 ITRF2014 给 出 的 当前 时 
刻 坐 标 进 行 比较 ， 最 后 计算 7 参数 ， 公 式 同 (2). 

基于 中 俄 合作 的 背景 ， 对 伪 距 用 户 兼 容 互 操作 提出 了 要 求 。 采 用 SLR 坐标 估计 法 对 
BDSC 与 PZ-90 之 间 的 转换 7 参数 进行 了 比较 。 

下 面 给 出 不 同 参考 框架 之 间 的 转换 公式 ， 以 BDCS 和 PZ-90 为 例 。 
PZ-90 与 ITRF 的 转换 公式 : 


jami 


日 公式 (2) 可 得 到 


XIrRF = X 


trans 


+ (1 + Xx! ) x Xt x Xpz—90 . (6) 


scale 


将 公式 (6) 代入 到 公式 (2) 中 ， 可 得 到 : 
Xpdes = X trans (1 Xscale) x Xiot x (XI F (1 Ez X! oale) x X'ot x Xpz_90) 


trans S (7) 
y n" 1" 
= A rans (1 XXscale) x X iot x Xpz—90 


日 此 可 得 到 BDCS 和 PZ-90 的 转换 参数 : 


EH 


X” 


trans 


2 
A scale = Xscale + A scale F A scale x Xscale ’ (8) 


/1/ / 
Arot x X ot 


rot 一 


= (1 N Xscale) X Xrot x X rans +E X¢trans ) 


3 BDCS 实 现 精度 监测 


3.1 ”试验 条 件 

本 文 试验 需求 的 监测 站 需要 具备 稳定 接收 北斗 三 号 B11/B3I 信号 ， 且 历 元 可 用 的 北斗 
三 号 卫星 至 少 4 颗 。 经 过 筛选 后 选用 15 个 台 站 数据 ， 分 别 为 境内 站 wuh2, urum; 境外 站 
tuva, zim2, ulab, arht, met3, savo, borl, pove, ganp, moiu, sthl, pots, sgoc。 表 1 为 各 站 历 
元 可 用 北斗 三 号 卫星 数量 统计 。 图 1 为 监测 站 站 址 分 布 图 。 


表 1 历 元 可 视 北 斗 三 号 卫星 数量 统计 表 


站 名 ”可 用 北三 卫星 数 | 站 名 可 用 北三 卫星 数 | 站 名 可 用 北三 卫星 数 
AHRT 7 POTS 6 TUVA 4 
BOR1 4 POVE 5 ULAB 5 
GANP 4 SAVO 5 URUM 6 
MET3 7 SGOC 6 WUH2 5 
MOIU 6 STHL 7 ZIM2 4 
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3.2 ”广播 星 历 定位 方法 结果 
采用 2019 年 1 月 到 12 月 全 球 监测 站 BDS 卫星 B1/B3I 双 频 伪 距 相位 30 s 采样 数据 及 
广播 星 历 ， 计 算 单 天 监测 站 坐标 。 统 计 了 全 年 的 坐标 估计 值 与 ITRF 坐标 互 差 ， 并 每 天 计算 


一 组 7 参数 。 


图 1 


监测 站 站 址 分 布 


采用 B11/B3I 双 频 组 合 ; 


因为 北斗 三 号 卫星 与 北斗 二 号 
与 北斗 二 号 、 三 号 联合 两 种 处 理 策 略 相 比较 。 第 一 种 全 部 使 有 
第 二 种 采用 北斗 二 号 /北斗 三 号 混合 星 
二 号 采用 B11/B3I 和 BI1I/B2I 双 频 组 合 ， 北 斗 三 号 为 B11/B3I 双 频 组 合 。 


星空 间 信号 精度 存在 差异 ， 所 以 我 们 仅 对 北斗 三 号 卫星 
北斗 三 号 了 卫星， 为 C19-C37， 


座 ， 为 C01-C37， 其 中 北斗 


北斗 三 号 处 理 各 站 2019 年 全 年 的 坐标 估计 值 与 TTRF 互 差 序列 的 统计 结果 ， 见 表 2。 
#2 ”北斗 三 号 广播 星 历 定位 坐标 互 差 


卫 方 向 /cm N 方向 /cm U 方向 /cm 

站 名 å mS oa JZ BE >Z DE JZ 
MET3 10503010 | —3.07 1.86 —17.91 1.14 4.66 3.61 
GANP 11515001 0.82 2.97 0.90 0.92 —7.31 3.52 
BOR1 12205002 1.28 3.08 0.65 1.38 一 7.75 4.69 
ZIM2 14001008 1.07 3.31 0.12 1.12 —9.00 5.29 
POTS 14106003 | —0.22 2.05 1.70 1.10 —14.45 3.99 
WUH2 21602007 | —0.06 1.34 1.14 0.48 一 2.22 2.48 
URUM 21612001 —0.28 1.39 1.35 0.31 2.26 2.41 
SGOC 23501003 | —4.30 2.02 0.10 0.35 一 5.00 3.45 
ULAB 24201001 一 1.23 1.43 1.60 0.49 一 2.80 2.18 
STHL 30606003 | —6.88 2.85 一 1.17 0.74 8.42 4.41 
MOIU 33204001 | —4.75 1.70 —0.39 0.40 4.87 2.60 
TUVA 51101001 一 1.61 32.61 一 1.25 3.13 5.44 16.91 
ARHT 66001005 3.12 1.53 一 2.31 1.02 4.68 5.54 
SAVO 93235000 | 一 4.98 2.75 0.42 0.44 —5.13 3.46 
POVE 93780000 | —3.70 3.60 0.43 0.50 一 4.92 4.12 
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我 们 利用 2019 年 2 月 北斗 三 号 数据 ， 通 过 广播 星 历 定位 法 扣除 卫星 钟 差 、 接 收 机 钟 差 、 


大 气 延 迟 以 及 相位 模糊 度 等 误差 ，3 
一 0.2 ~ 0.2m 之 间 。 


F 得 到 O 一 C 时 间 序 列 ， 结 果 如 图 2 Aras, 
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2019 年 年 积 日 


2 ”北斗 三 号 广播 星 历 定位 法 残 差 序列 


表 3、 图 3 给 出 了 北斗 三 号 广播 星 历 定位 法 的 转换 7 参数 结果 。 


表 3 北斗 三 号 广播 星 历 定位 法 7 参数 结果 


平均 残 差 在 


T,/em T,/em T./em R,/mas R,/mas Rs/mas 尺度 因子 /10-9 
参数 3.3 —0.7 3.9 一 2.19 0.08 一 0.43 0.76 
中 误差 2.0 5.6 4.5 1.88 1.52 1.54 0.53 
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图 3 


北斗 三 号 广播 星 历 定位 法 得 到 的 平移 参数 、 旋 转 参 数 和 尺度 因子 随时 间 的 变化 
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北斗 二 号 /三 号 联合 处 理 各 站 2019 年 全 年 的 坐标 值 与 ITRF 互 差 序列 见 表 4。 


表 4 北斗 二 号 /三 号 联合 广播 星 历 定位 坐标 互 差 
也 方向 /cm N 方向 /cm U 方向 /cm 
均值 ”方差 均值 方差 均值 方差 
MET3 10503010 | —1.47 0.94 —5.16 0.90 5.25 3.45 
GANP 11515001 | —1.28 2.06 8.90 6.96 —9.63 7.31 
BOR1 12205002 | —4.90 6.61 5.66 3.15 —8.02 2.22 
ZIM2 14001008 | —6.12 7.42 5.40 3.39 一 8.77 5.74 
POTS 14106003 0.21 0.85 2.66 0.63 一 9.27 5.56 
WUH2 21602007 2.14 0.97 2.72 0.81 一 4.53 7.64 
URUM 21612001 2.84 1.04 3.21 0.61 9.18 10.09 
SGOC 23501003 | —2.75 1.10 —1.49 0.79 —20.93 9.36 
ULAB 24201001 0.67 1.44 3.12 0.72 10.75 11.41 
STHL 30606003 1.20 2.92 0.43 0.44 5.86 6.73 
MOIU 33204001 | —1.47 1.02 —0.46 0.55 8.65 9.07 
TUVA 51101001 5.48 3.88 6.66 1.73 —16.51 7.38 
ARHT 66001005 5.13 3.18 —5.00 5.33 12.36 5.73 
SAVO 93235000 | —7.05 4.72 1.87 6.30 —8.85 8.96 
POVE 93780000 | —5.55 2.65 —2.43 3.66 一 7.24 5.36 


站 名 编号 


我 们 利用 2019 年 2 月 北斗 二 号 /三 号 数据 ， 通 过 广播 星 历 定位 法 扣除 卫星 钟 差 、 接 收 机 
钟 差 、 大 气 延 迟 以 及 相位 模糊 度 等 误差 ， 并 得 到 0 一 C 时 间 序 列 ， 结 果 如 图 4 所 示 ， 北 斗 三 
号 平均 残 差 在 -0.5 ~ 0.5 m 之 间 ， 北 斗 二 号 平均 残 差 在 -0.7 ~ 0.7 m 之 间 。 


ms 


2019 年 年 积 


2019 年 年 积 E 


图 4 北斗 二 号 /三 号 卫星 定位 法 北 二 /三 卫星 残 差 序列 


由 于 北斗 二 号 是 区 域 服务 系统 ， 没 有 境外 观测 和 注入 ， 因 此 进行 全 球 评估 时 采用 的 是 北 
斗 二 号 长 期 预报 产品 。 而 北斗 三 号 使 用 星 间 链 路 观测 ， 地 面 站 结合 星 地 星 间 数据 联合 定 轨 ， 
其 轨道 精度 较 高 。 因 此 相对 于 北斗 二 号 ， 北 斗 三 号 卫星 空间 信号 精度 较 高 ， 见 表 5 及 图 5。 
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表 5 北斗 二 号 /三 号 联合 广播 星 历 定位 法 7 参数 结果 
T,/om Ty/em T:/cm Rz/mas Ry/mas R.z/mas 尺度 因子 /1079 


参数 7.5 7.3 8.5 一 2.65 一 2.14 0.76 一 0.23 
中 误差 5.9 9.6 9.7 3.33 2.52 1.45 1.69 
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注 ，a) 平移 参数 ，b) 旋转 参数 ，c) 尺度 因子 。 
图 5 联合 广播 星 历 定 位 法 得 到 的 平移 人 参数、 旋转 参数 和 尺度 因子 随时 间 的 变化 


3.3” 星 历 比 对 法 结果 

采用 2019 年 1 月 3 日 到 1 月 27 日 总 计 25 d 的 精密 星 历 及 广播 星 历 ， 精 密 星 历 为 中 国 
科学 院 上 海天 文 台 事后 处 理 精密 星 历 。 将 广播 星 历 与 精密 星 历 相 拟 合 可 计算 出 卫星 在 BDCS 
框架 下 和 ITRF 框架 下 的 坐标 ， 并 统计 25 d 的 坐标 互 差 ， 每 天 计算 一 组 7 参数 。 
N 将 上 海天 文 台北 斗 事后 精密 星 历 与 广播 星 历 进行 比较 ， 计 算 单 天 7 参数， 并 计算 一 个 月 
~ 平均 值 ， 结 果 如 表 6 及 图 6 所 示 。 


表 6 BDS 广 播 星 历 与 精密 星 历 7 参数 结果 
T,/om Ty/em T:/em Rz/mas Ry/mas R.z/mas 尺度 因子 /10-9 


= 参数 一 3.4 0.5 3.5 0.01 0.40 0.04 —0.026 
kon 中 误差 1.6 0.5 0.9 1.10 0.81 0.08 0.04 


3.4 SLR 坐标 估算 法 结果 

采用 SLR 测量 数据 以 及 广播 星 历 ， 分 别 计 算 BDS 和 GLONASS 系统 的 SLR A ihih 
标 ， 结 果 见 表 7, 8。 对 于 BDS 结果 ，4 颗 北 斗 三 号 卫星 ， 共 3330 个 标准 点 数据 ; 对 于 
GLONASS 结果 ，21 颗 在 轨 卫 星 ， 共 4317 个 标准 点 数据 。 

我 们 采用 SLR 估算 法 解 得 2019 年 全 年 的 残 差 时 间 序 列 ( 见 图 7)， 表 明 该 方法 下 BDS 
的 空间 信和 号 精度 的 残 差 在 -0.15 ~ 0.15 m 范围 内 。 

采用 SLR 测量 数据 以 及 广播 星 历 ， 计 算 SLR 台 站 坐标 。 数 据 为 2019 年 BDS 全 年 
ILRS 激光 联 测 数据 ， 由 于 SLR 数据 采样 较 少 ， 一 个 月 估计 一 组 台 站 坐标 ， 计 算 了 BDCS 
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iE: a) 平移 参数 ，b) 旋转 参数 ，c) 尺度 因子 。 
6 ” 星 历 比 对 法 得 到 的 平移 参数 、 旋 转 参 数 和 尺度 因子 随时 间 的 变化 


表 7 GLONASS SLR ,坐标 估计 互 差 


编号 EE 方向 /cm NN 方向 /cm U 方 向 /cm 

均值 方差 均值 方差 均值 ”方差 

18685901 5.27 1.30 —5.97 1.24 —13.72 1.07 
78393402 一 3.81 2.05 一 0.12 2.25 一 9.23 1.58 
78403501 5.44 1.71 —0.88 1.63 —8.35 1.39 
78418701 一 2.55 4.67 0.31 3.86 一 9.41 2.64 
72496102 11.71 2.80 一 2.53 2.29 一 14.27 2.06 
78212801 2.86 2.25 —9.07 1.88 —3.84 1.61 
18799401 —0.63 2.56 一 1.52 3.14 一 5.92 2.07 
18734901 一 11.83 6.50 一 6.70 6.77 一 7.01 4.75 
18879701 3.54 10.45 一 7.85 16.06 一 8.13 11.8 
= 70900513 3.49 1.22 —10.55 1.31 一 1.14 0.92 
18900901 一 0.16 5.38 一 16.03 4.40 一 6.23 4.51 
78259001 —10.64 2.69 —2.08 3.09 18.34 2.10 
75010602 257.90 5.55 —152.46 5.71 —94.33 7.95 
71240802 | —256.41 5.55 —152.46 5.71 —94.33 7.95 
79417701 11.00 2.65 3.65 3.79 —12.62 2.19 
78457801 —13.33 2.62 —1.39 2.34 一 5.42 2.22 
71050725 58.13 3.28 一 228.27 3.43 187.27 3.18 
71100412 17.98 6.28 4.52 3.94 一 15.29 3.67 
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表 8 2019 年 BDS SLR 坐标 估计 互 差 
编号 了 方向 /cm N 方 向 /cm U 方 向 /cm 
均值 FH 均值 Fe 均值 方差 
70900513 | -6.73 5.09 -2.20 8.78 —3.49 7.09 
78457801 | 4.64 5.40 -1.23 3.71 3.07 6.76 
78393402 | 1.90 15.17 0.32 12.86 —8.45 15.12 
78403501 | 2.33 643 -0.33 3.33 2.06 8.03 
88341001 | 1.77 4.30 437 7.78 6.51 5.42 
71100412 | -4.18 19.95 6.89 18.30 3.73 17.69 
72496102 | 4.22 9.45 11.50 15.74 -3.37 16.62 


7 SLR 坐标 估算 法 残 差 序列 


与 ITRF 之 间 的 转换 7 参数。 此 外 ， 采 用 GLONASS 2019 年 1 月 ILRS 激光 联 测 数据 ， 计 
@ 7 PZ-90 与 ITRF， 以 及 BDCS 与 PZ-90 之 间 的 转换 7 参数， 结果 见 表 9 及 图 8。 


表 9 SLR 坐标 估算 法 7 参数 结果 


Ts/cm Ty/em T/cm Rz/mas Ry/mas Rs/mas 尺度 因子 /10- 
参数 一 4.8 3.4 2.2 2.46 一 2.60 0.78 0.004 
BDS a 
中 误差 9.6 4.3 9.8 4.52 2.99 4.72 0.79 
参数 | —0.045 1.6 -1.2 -1.31 一 0.99 一 0.94 一 1.802 
GIONASS a 
中 误差 2.1 2.1 2.0 0.78 0.71 0.87 0.3 
参数 一 0.3 17 一 6.5 一 4.51 0.01 —0.85 —2.786 
BDS-GLONASS K 
中 误差 5.7 5.4 5.2 1.91 2.05 2.24 0.8 
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VE: a) 平移 参数 ，b) 旋转 参数 ，c) 尺度 因子 。 


图 8 SLR 坐标 估算 法 得 到 的 平移 参数 、 旋 转 参 数 和 尺度 因子 随时 间 的 变化 


3.5 ”结果 分 析 

本 文采 用 广播 星 历 定位 法 、 广 播 星 历 与 精密 星 历 比较 以 及 SLR 台 站 坐标 估算 三 种 方法 
对 BDCS 实现 精度 进行 计算 。 第 一 种 方法 分 别 用 北斗 二 号 、 北 斗 三 号 混合 星座 和 北斗 三 号 
进行 精密 坐标 估算 。 北 斗 三 号 7 参数 转换 结果 平移 参数 为 3.3 cm, —0.7 cm, 3.9 cm; Wel 
数 为 -2.19 mas, 0.08 mas, 一 0.43 mas; 尺度 因子 为 7.6 x 10-83。 北 斗 二 号 和 北斗 三 号 混合 
星座 7 参数 转换 结果 平移 参数 为 7.5 cm, 7.3 cm, 8.5 cm; 旋转 参数 为 -2.65 mas, 一 2.14 mas, 
0.76 mas; 尺度 因子 为 ~2.3 x 10-10。 

第 二 种 方法 将 广播 星 历 与 精密 星 历 进行 比较 。 转 换 结果 平移 参数 为 -3.4 cm, 0.5 cm, 
3.5 cm; 旋转 参数 为 0.01 mas, 0.4 mas, 0.04 mas; 尺度 因子 为 -2.6 x 10-11。 

第 三 种 方法 采用 精密 星 历 和 SLR 数据 计算 台 站 坐标 。BDS 转换 结果 平移 参数 为 
一 4.8 cm, 3.4 cm, 2.2 cm; 旋转 参数 为 2.5 mas, —2.6 mas, 0.8 mas; 尺度 因子 为 4 x 10-12。 
GLONASS 转换 结果 平移 参数 为 -0.045 cm, 1.6 cm, -1.2 cm; 旋转 参数 为 -1.31 mas, 
一 0.99 mas, 一 0.94 mas; 尺度 因子 为 —1.802 x 107°. 

BDS 与 GCLONASS 转换 结果 平移 参数 为 -0.3 cm, 1.7 cm, —6.5 cm; 旋转 参数 为 
一 4.51 mas, 0.01 mas, —0.85 mas; 尺度 因子 为 —2.786 x 107°. 
通过 上 述 计 算 结 果 可 以 看 出 ， 三 种 方法 计算 的 了 参数 的 形式 误差 均 在 厘米 量 级 ， 因 此 
BDCS 与 ITRF 对 齐 精度 在 厘米 量 级 。 

对 比 上 述 三 种 方法 ， 由 于 星 历 比 对 法 和 SLR 坐标 估算 法 引入 了 事后 精密 产品 ， 监 测 精 
度 比 广播 星 历 定位 法 高 。 北 斗 卫 星 广播 星 历 包含 轨道 误差 、PCO 误差 等 ， 因 此 所 得 精度 不 
高 ， 会 对 参考 框架 实现 精度 监测 造成 一 定 的 影响 ， 但 随 着 广播 星 历 精度 的 提高 ， 其 监测 精度 
也 会 逐渐 提高 。 


ChinaXiv 合 作 期 刊 


116 天 文学 进展 39 卷 


4 总 结 


本 文 分 别 给 出 了 三 种 方法 来 对 BDCS 的 实现 精度 进行 计算 ， 并 采用 了 2019 年 全 年 的 数 
据 进 行 试验 ， 结 论 如 下 。 

(1) 通过 上 述 三 种 方法 对 实测 数据 进行 计算 ， 均 可 满足 RNSS 米 级 用 户 兼 容 互 操作 的 需 
求 。 广 播 星 历 定位 法 可 以 进行 参考 框架 的 实时 监测 ， 但 广播 星 历 精度 决定 了 监测 精度 ， 轨道 
比 对 法 通过 引入 事后 精密 产品 来 提高 监测 精度 ， 但 时 效 性 不 强 ，SLR, 坐标 估算 法 所 得 的 精 
度 是 三 种 方法 中 最 高 的 ， 但 目前 激光 联 测 数据 量 较 少 ， 因 此 只 能 进行 长 期 的 监测 结果 分 析 。 

(2) 本 文 所 选取 的 IGS 站 均 需 能 够 接收 到 北斗 三 号 卫星 BII/B3I 频 点 ， 因 此 未 选取 太 多 
站 参与 本 次 试验 。 北 斗 二 号 和 北斗 三 号 联合 处 理 的 方法 与 仅 北 斗 三 号 处 理 的 结果 在 厘米 量 
级 一 致 ， 因 此 在 后 续 的 工作 中 ， 可 以 加 入 更 多 以 IGS 核心 站 为 主 的 地 面 监测 站 进行 全 球 网 
解 计算 ， 以 进行 参考 框架 监测 ， 有 助 于 监测 结果 的 可 靠 性 。 
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Abstract: From the perspective of RNSS meter-level services, BDCS implementation ac- 
curacy evaluation methods are given, including broadcast ephemeris positioning method, 
precise ephemeris and broadcast ephemeris comparison method, and satellite laser ranging 
(Satellite Laser Ranging, SLR) coordinate estimation method. The broadcast ephemeris 
positioning method is suitable for users to realize real-time monitoring of the accuracy of 
the reference frame, but the monitoring accuracy is limited by the accuracy of the broadcast 
ephemeris. The orbital comparison method introduces post-precision products with high 
monitoring accuracy, but its timeliness is relatively poor. The SLR coordinate estimation 
method is less affected by the ionospheric delay and is not affected by satellite clock errors. 
It has the highest accuracy among the three methods. However, due to the small amount of 
laser data, it cannot be used for short-term monitoring. The results show that the accuracy 
of the three methods to evaluate the alignment of BDCS and ITRF is at the centimeter 
level. Finally, in the context of Sino-Russian satellite navigation system cooperation, the 
alignment accuracy of BDCS and PZ90 was evaluated. The results show that the two are 
aligned on the centimeter level, which can meet the application requirements of meter-level 
positioning users. With the improvement of Beidou-3 broadcast ephemeris accuracy and the 
development of SLR, the evaluation accuracy of the above three methods will be further 


improved. 
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